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Abstract 

We show that a straightforward rewrite of a known minimal polynomial algorithm 
yields a simpler version of a recent algorithm of A. Salagean. 
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1 Introduction 

Let K be a field, n > 1 and s = (si, . . . , s n ) be a finite sequence over K. The Berlekamp- 
Massey (BM) algorithm computes an LFSR of shortest length L and a feedback polyno- 
mial F G K[x] generating s, vacuously if L = n pp. 

We begin with the approach and basics of [2J. Multiplication makes Laurent series 
in x~ l into a .fT [a;] -module and power series with non-zero annihilator ideal correspond 
to linear recurring sequences: they have a non-zero 'characteristic polynomial' (c.p.) [21 
Section 2]. For finite sequences, we work with Laurent polynomials and C G K[x] is a 
c.p. of s if 

Co ■ s 3 . d + ■ ■ ■ + C d ■ Sj = (1) 

for d+ 1 < j < n, where d = deg(C) > [21 Definition 2.7, Proposition 2.8]. Any C with 
d > n is vacuously a c.p. of s. 

A c.p. C of s is a minimal polynomial of s if deg(C) = min{deg(D) : D is a c.p. of s} 
and the linear complexity of s is the degree of any minimal polynomial of s [21 Definition 
3.1], [31 Definition 2.2]. For example, D = x L ~ deg( - pS) F* is a c.p. of s — as usual, F* is the 
reciprocal of F — and D is a minimal polynomial of s since D* = F and deg(D) = L. 
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As far as we know, Algorithm 4.2 of [2J (Algorithm 3.1 below) was the first algorithm 
to compute a minimal polynomial of s iteratively. In fact, it is valid for finite sequences 
over a commutative unital integral domain. 

Algorithm 2.2 of [4J also computes a minimal polynomial of s. We show that these two 
algorithms are closely related: a straightforward rewrite of the former using the notation 
of [I] yields the latter, except that we initialise a polynomial to instead of 1. Further, 
the rewrite uses fewer variables and is simpler. See also Remark 3.2 (iv). 

We note that [2J and [3] (an expository version of [2j) were referred to in [51 Introduc- 
tion] . 

2 The Inductive Construction 

The whole process of Algorithm 3.1 is best explained in terms of the inductive construction 
of a minimal polynomial of s which was derived from first principles in [2] . 

2.1 The Naive Version 

A natural choice for is 1 if si = and x otherwise; is certainly a c.p. of minimal 
degree. Now assume inductively that 2 < i < n and we have c.p.'s C^> for (si, . . . ,Sj) 
where 1 < j < i — 1. From Equation ([1]), C^~^ is a c.p. of (si, . . . , Si) if and only if the 

discrepancy 

c i-i = d(C^) = " 'i-H-*-! e * 

5=0 

is zero, where d^ x = deg{C^) Definition 2.10], cf. [1, Equation (10)]. If q_i = 
then clearly c? i _ 1 is minimal. But if Cj_i 7^ 0, a new c.p. is needed. We use an index a;_i 
suggested by [H Equation (11)] 

a = aj_i = max {?' : gL < cL 1} (2) 

l<7<«-2 

to index a previous c.p. [21 Definition 3.12]. The exponent e^-i = 2dj_i — i and are 
from [21 Proposition 4.1]: 

r c a • C^- 1 ) - q_i • x +e C( a ) if e = ei _! > 
C (i) = ^ (3) 
[ c a • x _e C^ -1 ^ — Cj_i ■ otherwise. 

Now a = aj_i and Equation ([3]) require i > 3 and > di. For 2 > 2 and di-i = d%, we 
complete our construction by 

cW = \ x l if si = (4) 

1 sx • 2 C^ ^ — Cj_i otherwise. 
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The case si = occurs when s has i — 1 > 1 leading zeroes. 

Part (i) of the following proposition is an analogue of [TJ Equation (15)]. For Part (ii), 
see [3], Propositions 4.3, 4.5]. 

Proposition 2.1 For 2 < % < n, if q_i 7^ then (i) di = max{dj_i,i — and (ii) 
is a c.p. of (si, . . . , Si). 

An analogue of [TJ Lemma 1] establishes minimality. 

Lemma 2.2 ([21 Theorems 3.8, 3.13] or [31 Lemma 5.2]) Let 2 < z < n, / be a c.p. of 
(si, . . . , Sj_i) and d(/) 7^ 0. If # is a c.p. of (si, . . . , S,), then deg(g) > i - deg(/). 

We now have a naive inductive construction for a minimal polynomial of s and illus- 
trate it with two binary examples. 

Example 2.3 Consider the subsequence s = (0, 1, 1, 0) of [H Table I]. We have = 1 
as si = and C\ = s 2 7^ 0. As i < 3, we apply Equation (J3J): = x 2 (there is one 
leading zero). Now i = 3, c 2 = s 3 7^ and <i 2 = 2 > = 0, so that a 2 = 1. Equation (j3J) 
applies with exponent e 2 = 2- 2 — 3 = 1 giving = + x 1 C^ = x 2 + x. Finally, 
c 3 = s 4 + s 3 7^ and d 3 = 2 > d\ = 0, so that a 3 = a 2 = 1. Equation (j3~J) applies again 
with exponent e 3 = 2 ■ 2 - 4 = and C (4) = C (3) + a;°C (1) = x 2 + x + 1. 

Example 2.4 Let s = (1,1,0,0). Clearly CM = x and ci = s 2 7^ 0, so C< 2 > = 
x°C^ + 1 = x + 1 from Equation (j3J). Next, c 2 = s 3 + s 2 7^ and Equation (j3J obtains 
as d 2 = di = 1, giving = a; 1 ^ 2 ) + c 2 = x 2 + a; + 1. For i = 4, c 3 = s 4 + s 3 + s 2 7^ 0, 
ofa = 2 > g?i = 1 and a 3 = 2. Equation ([3]) with exponent e 3 = 2- 2 — 4 = gives 
=C®+x°C<® =x 2 . 



2.2 The Refined Version 

The naive construction can be refined in three ways. Firstly, by noting that if ej_i > 
then a, = aj_i since d^\ > i — i.e. di = di-i by Proposition 2.1. But if e^_i < 0, 
di = i — di-\ > c?j_i and dj = i — 1. 

Secondly, we update ej_i = 2dj_i — i (and so avoid using c?j_i) as follows: it is trivial 
that e, = ej_i — 1 if q_i = and easy to check that e, = | 1 1 — 1 if Cj_i 7^ 0. 

Thirdly, we change the inductive basis by introducing artificial values C^ ' = 1 (which 
is only a c.p. if si = 0), C^ ) = (which is not a c.p. by definition) and c ao = 1 (which is 
not a discrepancy). For definiteness, we take = —1. Then remarkably (i) eo = — 1 (ii) 
Equation (j3~J) accomodates all three cases and (iii) the updating of and ej_i remains 
valid. We state this formally. 

Theorem 2.5 (0) Put a = e = -1, = 0, c_i = 1 and = 1. For 1 < i < 71, 

define C^, a,, ej by 

(i) if Cj_i = 0, = C^" 1 ), aj = aj_i and e, = ej_i — 1 
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(ii) if Ci-i 7^ let CW be as in Equation (EJ), = aj_i if ej_i > 0, a, = z — 1 if ej_i < 
and ti = |ej_i| — 1. 

Then for 1 < i < n, is a minimal polynomial of (si, . . . , Sj) and deg(CW) = £i+ ' + . 

The BM algorithm decodes Reed-Solomon, Goppa and negacyclic codes [6], [7] and 
has been extended to multiple sequences [8]. For similar applications and extensions of 
Theorem 2.5, see [3j Section 8], [9]. 

3 The Algorithm 

We deduce Algorithm 3.1 from Theorem 2.5. Firstly, we can dispense with the indices a\ 
if we define B^ = C^ ai ' and scalars bi = d(C^) for < i < n — 1 provided we update 
B^~^ and when e;_i < 0. Secondly, we replace by e '~ 2 1+t in Finally, since 
only current values are used we can suppress all indices — provided we keep a copy T of 
C^-V to update when e^i < 0. 

Algorithm 3.1 ([21 Algorithm 4.2] rewritten) 

Input: n > 1 and a sequence s = (si, . . . , s n ) over a field fT. 

Output: C, a monic minimal polynomial for s. 

begin £? 0; b <- 1; C <- 1; e ^ 1; 

for i = 1 to n do 

if c ^ then if e > then C <- b ■ C - c • x e £; 

else T <- C; e < e; 

C <- 6 • x e C - c ■ B; 

B <-T;b<-c; endif 
endif 

e <— e — 1; 
endfor 

return (C/b) end. 

Remark 3.2 We obtain [H Algorithm 2.2] except that 

(i) we could (but do not) make each C monic 

(ii) in the notation of [4, p. 4696], we do not keep track of 'm' or deg(S) to recompute 
v = N — m — (deg(C) — deg(B)) at each iteration (where < N < n — 1) 

(iii) in fact, (a) N corresponds to % — 1 (b) if there are N leading zeroes v = N +1 (c) 
if deg(C) > deg(B) then deg(C) = m + l-deg(5) and v = N-m- (deg(C) -deg(B)) = 
N + 1 — 2 deg(C); thus v corresponds to — e 

IV = and some initial minimal polynomials will differ; with B^ — ^ = 1, 
= x 1 — Cj_i when there are % — 1 > leading zeroes, Theorem 2.5 remains valid, the 
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algorithms are equivalent and their outputs are identical. 



Remarks 3.3 (added in proof) 

(i) Theorem 2.4 of |4J on the set of minimal polynomials of s is an immediate con- 
sequence of pi Theorem 4.16]: simply replace d in Theorem 4.16, loc. cit. by — v = 
deg(C) — deg(_B) — (n — m) and take first components. 

(ii) In [10], we calculate the monic reciprocals of the minimal polynomials of The- 
orem 2.5. This readily yields an algorithm similar to the BM algorithm, except that 
we do not calculate the 'lengths' L,; we update the exponents as above and obtain L n as 
£a±p+I . This article also contains proofs for the results of Section [2] using similar notation. 

The author thanks the anonymous reviewers for their helpful comments and sugges- 
tions. 
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